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Abstract 


Some  results  of  high  speed  interferometric  measurements  on  dynamically 
propagating  inteifacial  cracks  are  presented  PMMA/Steel  bimaterial  specimens  of  the 
one  point  bend  type  are  used  They  are  inq)act  loaded  using  eitho*  a  drop  weight  tower 
device  or  a  high  speed  gas  gun.  This  results  in  two  distinct  groups  of  experimmts  at 
different  loading  rates.  In  all  cases  very  hig^  crack  propagation  speeds  are  seen.  In  gas 
gun  experiments,  terminal  crack  tip  speeds  of  up  to  1.5Cs^^®^,  where  ^  jjjg 

shear  wave  speed  of  PMMA,  are  measured  In  addition,  high  accelerations  (»10^g,  where 
g  is  die  acceleration  of  gravity)  are  observed  and  reported.  Theoretically  predicted  near 
tip  fields  are  used  to  extract  eiqietimental  values  of  the  dynamic  con^lex  stress  intensity 
factor  histories  in  each  test  Using  die  fitted  histories  of  the  complex  stress  intensity 
factor,  a  dynamic  crack  growth  criterion  relating  the  energy  release  rate  to  phase  angle 
and  crack  t^  velocity  is  proposed  and  discussed 
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1  Introduction 

Tlw  increased  use  of  con^osite  materials  in  engineering  applications  in  recent 
decades  has  become  an  incentive  for  extensive  research,  both  basic  and  applied,  into  the 
various  failure  rtrades  of  such  materials.  It  was  soon  recognized  that  the  macroscopic 
perfomumce  of  composite  materials  was  intimately  related  to  effects  occurring  at  the 
interface  betwem  the  different  phases  of  the  composite.  This  led  to  a  recent  explosion  of 
research  in  interfacial  firacture  effects  in  bimaterial  systems. 

The  first  study  of  interfacial  firacture  was  made  by  Williams  (1959).  He  found  an 
oscillatory  stress  singularity  present  at  the  tip  of  a  crack  lying  on  an  interface  between  two 
elastic  solids.  Until  recently,  only  a  few  researchers  had  provided  some  explicit  solutions 
of  the  stress  field  in  particular  interfacial  fracture  problems  (Sih  and  Rice  (1964),  England 
(1965),  Erdogan  (1965),  and  Rice  and  Sih  (1965)).  In  all  these  solutions,  asymptotically 
an  oscillatory  stress  singularity  was  found,  confirming  Williams  (1959)  result  The  recent 
resurgence  of  work  on  interface  firacture  has  resulted  in  many  more  contributions.  Most 
notable  are  those  of  Rice  (1988),  and  Shih  (1991).  Up  to  this  point  most  of  the  research 
activity  on  interface  fiacture  is  on  issues  regarding  the  quasi-static  failure  of  bimaterial 
systems.  To  folly  understand  the  catastrophic  nature  of  composite  fracture,  which  is 
important  in  many  engineering  applications,  a  good  understanding  of  dyruanic  interfacial 
fiacture  mechanics  is  also  necessary.  Due  to  the  conq>lexity  of  this  problem,  however,  only 
a  few  theoretical  investigations  have  been  made  so  far  (Goldshtein  (1967),  Brock  and 
Achenbach  (1973),  Wilis  (1971,  1973),  Atkinson  (1977),  Yang  et  al.  (1991)  and  Deng 
(1993)).  To  our  knowledge,  only  two  preliminary  experimental  studies  on  this  subject  are 
currently  availaUe  (T^ur  and  Rosakis  (1991)  and  liu  et  al.  (1993)). 

In  tile  current  work  we  investigate  the,  essentially  unexplored,  area  of  experimental 
ctynamic  interfacial  firacture  mechanics.  The  general  objective  sought  by  tiiis  investigation  is 
the  undorstanding  of  the  highly  transient  dynamic  interfacial  crack  initiation,  growth  and 
arrest  in  bimatetial  systems.  In  particular,  in  the  context  of  the  present  paper,  it  is  intended 
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to  provide  a  sommaiy  of  recent  e}q)erimental  observations  on  the  issue  of  dynamic  interface 
fracture,  culminating  at  the  proposal  of  a  physically  based  and  accurate  criterion  for  crack 
initiation  and  growth  in  bimaterial  systems.  Use  of  high  speed  photography  and 
interferometric  optical  diagnostics  can  provide  real  time  measurements  of  mechanical  fields 
at  the  vicinity  of  initiating  and  propagating  interfacial  crack  dps.  The  particular  optical 
method  used  here  is  die  newly  developed  method  of  Coherent  Gradient  Sensing  (GGS). 
(Tippur  et  al.  (1991),  Rosakis  (1993))  Hiis  is  a  full  field  lateral  shearing  interferometer 
which  has  shown  great  promise  in  teal  time  imaging  of  dynamic  crack  tip  fields  in 
hcmiogeneous  materials. 


2  Experimental  procedure 

2.1  Coherent  Gradient  Sensing  interferometer  (CGS) 

The  method  of  CGS  was  first  proposed  by  Uppur  et  al.  (1991)  and  is  described  in 
detail  by  Rosakis  (1993).  The  set-up  of  the  mediod  in  a  transmission  configuration  for  the 
case  of  a  dynamically  loaded  bimaterial  specimen  is  shown  schematically  in  figure  1.  A 
coherent,  monochromatic,  collimated  laser  beam  is  incident  on  the  deforming  specimen. 
After  transmission  through  the  specimen  it  acquires  an  optical  path  diffnence  and  looses 
collimation.  The  resulting,  now  non-collimated,  beam  passes  through  two  line  diffraction 
gratings  Gi  and  G2  of  fine  pitch  p  (typically  40  lines/mm).  They  are  situated  a  distance  A 
apart  and  perform  a  shearing  of  the  incident  wave  front  as  described  below.  The  gratings’ 
output  intensity  is  transmitted  through  a  filtering  lens  L.  A  dii^ction  spot  pattern  is 
obtained  on  the  filtering  plane,  which  is  located  at  the  back  focal  plane  of  tiie  lens  L.  At  this 
plane  all  but  one  difiraction  orders  are  blocked.  The  one  remaining  diffraction  spot  (either 
of  ±1  orders),  shown  in  figure  1  as  the  open  circle  on  the  filtering  plane,  is  imaged  to 


-3- 


pcoduce  an  interfeitnice  pattern.  For  the  case  of  a  dynamic  expaiment  die  imaging  medium 
is  high  qieed  camera  focused  on  the  specimen. 

If  S(xiJC2)  is  the  optical  path  change  introduced  by  the  deforming  specimen  in  the 
transmitted  laser  beam,  then  it  has  been  shown  by  Tippur  et  al.  (1991)  diat  constructive 
inletfnence  will  occur  on  the  image  plane  when 


d(S(x,,X2))  _  mp 
dXj  A  ’ 


m  =  0,±/,±2,...  . 


(1) 


If  die  grating  lines  are  parallel  to  die  xi  direction  then  it  can  be  shown  that  die  condition  for 
ccHistmctive  interference  becomes. 


^S(xj,X2))  _  /ip 

ax2  ”4’ 


n  =  0,±/,±2,...  . 


(2) 


For  solid  mechanics  qiplications,  it  is  desirable  to  relate  the  quantity  S(X]^2)  to 
the  stress  state  in  the  deforming  specimra.  Such  a  relation  is  discussed  in  detail  in  Rosakis 
(1993)  and  will  not  be  presented  here.  For  a  linearly  elastic,  isotropic  material  under  plane 
stress  conditions,  see  Lee  and  Rosakis  (1993)  fOT  the  extent  of  the  plane  stress  region  in  a 
cracked  bimaterial  specimen,  5(X7,X2)  can  be  expressed  as, 

S(X„X2)  «  C^O„{X2.X2)-^ ^22{Xt.X2)],  (3) 

where  Ccr  is  a  stress  optical  coefficient  for  the  material,  /i  is  the  specimen  thickness  and 
and  &22  thickness  averages  of  the  in-plane  stress  con^nents  in  the  plate. 

Thus  for  points  outside  the  near  tip  three  dimensional  region  the  CGS  patterns 
assume  a  sinqile  interpretation  in  terms  of  two  dimensional  stress  field  approximations.  In 
particular,  equations  (1)  and  (2)  in  conjunction  with  equation  (3)  now  indicate  that  the 


-4- 


fringfis  obtained  from  regions  surrounding  the  duee  dimensional  zone  can  be  related  to  the 
in>plane  gradients  of  aj,  +  <722  as  follows 


=  =  m.n=0.±I.±l . 


A’  dbc,  A 


(4) 


^tdiere  in  die  case  of  transmission  Ca  is  the  stress  qptical  coefficient  of  FMMA. 


IJl  Binuterial  Experiments 

To  be  able  to  use  a  transmission  CXjS  arrangement  in  our  acperimentadon,  one  side 
of  the  bond,  the  more  compliant  side,  was  chosoi  to  be  Plexi^as  (Poly-Mediylmetfaacrylate 
or  PMMA).  This  is  a  transparent  polymer  that  is  easily  machined  and  handled.  At  room 
temperature  it  behaves  in  a  britde  fashion.  The  other  side  of  the  bond  was  chosen  as  eidier 
6061-T6  aluminum  or  AISI 4340  steel,  both  of  which  ate  considerably  stiffer  and  tougher 
than  Plexiglas.  Throughout  this  study,  the  PMMA  side  of  each  specimen  will  be  referred  to 
as  mateiial-l  and  the  metal  side  as  matetial-2.  So  with  reference  to  figure  2,  the  more 
coiiqiliant  material  occupies  the  positive  X2  plane.  The  mechanical  properties  of  the 
constituents  are  shown  in  table  1.  Table  1  also  shows  the  value  of  the  plane  stress  quasi¬ 
static  oscillatory  index,  defined  later  on  by  equation  (7),  for  the  bimaterial  combinations  of 
PMMA/aluminum  and  PMMA/stecl  Also  shown  is  the  value  of  for 
Plexiglas  bonded  to  a  rigid  substrate.  As  can  be  seen,  these  three  values  are  extremely 
close.  This  would  suggest  that  we  can  expect  intense  interfacial  effects  during  the  dynamic 
fracture  event  It  also  suggests  that  the  approach  of  some  researchers  (Hutchinson  (1989), 
Nakamura  (1991))  to  set  &=Q  invoking  small  material  mismatch  cannot  be  applied  in  our 
case. 

The  specimen  preparation  procedure  follows  the  lines  of  Tippur  and  Rosakis  (1991) 
and  Liu  et  aL  (1993)  and  will  not  be  described  here.  The  crucial  issue  of  bond  strength  has 
been  addressed  by  the  performance  of  botKi  calilxation  tests  by  Tippur  and  Rosakis  (1991). 
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These  tests  essentially  consist  of  static  tests  to  measure  the  fracture  toughness  of  specimens 
made  up  of  two  halves  of  PMMA  bonded  together  using  the  procedure  described  above.  It 
was  shown  that  the  firacture  toughness  of  the  bonded  PMMA  specimens  was  over  95%  of 
that  of  homogeneous  PMMA. 

The  biinaterial  specimens  used  in  the  dynamic  experiments  have  either  a  pre-cut 
edge  notch  or  a  shaip  pre-crack  of  length  2Smm  along  the  interface.  After  the  initiation 
event,  a  crack  propagates  dynamically  along  the  interface.  The  bimaterial  specimens  are 
impact  loaded  in  a  drop  wdght  tower  (Dynatup-8100A)  or  a  high  speed  gas  gun.  The 
transmission  GGS  technique  in  conjunction  widi  high  q>eed  photogr^hy  is  used  to  record 
dynamic  crack  tip  fields  around  the  notch  tip  (only  on  the  PMMA  side  of  course).  A 
rotating  muror  high  speed  camoa  (Gordin  model  330a)  is  used.  A  Spectra-Physics  Argon- 
ion  pulse  laser  (model  166)  is  used  as  the  light  source.  By  using  short  pulses  of  30ns 
duration  we  are  able  to  freeze  even  the  fastest  of  running  cracks  and  thus  produce  a  sharp 
interference  pattern  during  crack  growdi.  The  interframe  time  (controlled  by  the  interval 
between  pulses)  is  typically  set  at  Ips  for  a  total  recording  time  of  80^s.  The  laser  pulsing 
is  triggered  by  a  strain  gauge  on  the  specimen  that  smses  the  impact  of  the  drop  weight 
When  it  reaches  a  certain  threshold,  the  laser  starts  pulsing. 

True  symmetric  one  or  three  point  bend  loading  carmot  be  achieved  since  it  is 
extremely  difficult  to  tqrply  the  inq>act  load  exactly  on  the  interface,  which  is  very  thin.  In 
addition  since  the  wave  speeds  of  PMMA  and  steel  are  vastly  different  the  loading  history 
at  the  crack  ti^  would  be  con^ktely  different  if  die  specimen  were  impacted  on  the  PMMA 
or  die  steel  side.  Thus  it  was  chosen  to  impact  die  specimen  a  small  distance  (7mm)  into  die 
steel  side  of  the  bond.  A  sequence  of  high  speed  interferograms  from  a  PMMA/aluminum 
test  is  shown  in  figure  3  (from  Uppur  and  Rosalds  (1991)).  This  is  a  three  point  bend  test 
conducted  in  a  drop  weight  tower.  The  impact  speed  used  was  2  m/s.  When  the  crack 
initiates  (t=0ps),  intense  stress  waves  emanate  from  the  crack  tip.  These  waves  are  visible 
in  figure  3  as  discrete  kinks  in  othmwise  smooth  fringes  and  as  circular  lines  centered  at 


-6- 


points  along  dw  crack  Une.  This  observation  is  a  reliable  sign  of  a  highly  dynamic  event 
Indeed,  as  will  be  seen  later,  terminal  velocities  seem  to  be  about  90%  of  the  Rayleigh 
wave  speed  of  PMMA,  and  in  fact  in  some  tests  they  have  exceeded  it  Similar  waves  are 
visible  in  the  interferograms  of  figure  4.  These  are  fiom  a  PMMA/steel  drop  weight  tower 
one  bend  nperiment  at  an  impact  velocity  of  ^n/s. 


3.  Subsonic  and  intersonic  crack  growth  alcmg  bimaterial  interfaces 

3.1  Experimental  observations 

Drcp  weight  tower  experiments 

Figure  3  shows  a  sequence  of  dynamic  interferograms  from  a  PMMA/aluminum 
three  point  bend  test  (from  Tlppur  and  Rosaids  (1991)).  Time  t^Ops  corresponds  to  the 
time  of  crack  initiation.  Inteiferograms  at  negative  times  have  been  taken  before  initiation. 
In  this  figure  it  can  be  seen  that  between  t=-42ps  and  t=0iJs  the  fringe  pattern  surrounding 
the  initial  notch  tip  grows  in  size  and  also  rotates  with  req>ect  to  the  notch  line.  The  size  of 
the  fringe  pattern  at  each  time  dq>ends  on  the  magnitude  of  the  stress  field  around  the  tip. 
The  particular  orientation  of  each  lobe  of  the  pattern  is  dq)endent  on  the  relative  amounts  of 
opening  and  shearing  stresses  around  the  tip  (i.e.,  the  stress  mixity).  It  is  clearly  visible 
that  before  initiation  the  pattern  surrounding  the  notch  tip  changes  in  both  size  and 
orientation.  In  contrast,  after  crack  initiation  the  fringe  pattern  stays  relatively  uniform, 
signifying  that  some  fundamental  physical  quantities,  such  as  stresses  or  crack  face 
displacements,  must  remain  constant  during  crack  propagation.  Such  behavior  is  typical  of 
all  tests  performed  in  under  the  drop  weight  towo-  device. 

After  the  specimra  is  impacted,  on  the  metal  side,  a  compressive  stress  wave  travels 
the  vertical  width  of  the  specimen.  It  is  reflected  as  a  tension  wave  from  the  traction  free 
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end  on  the  bottom  of  the  plate.  Aldiough  impact  is  always  on  the  side  of  the  metal,  impact 
waves  are  also  generated  in  the  PMMA  side  because  of  the  existence  of  the  bond. 
However,  waves  in  the  metal  side  travel  much  faster  than  in  the  PMMA  side  so  the  wave 
generated  upon  inq>act  reaches  the  notch  tq)  from  the  metal  side  first  The  tip  is  thus  loaded 
by  subsequent  arrivals  of  reflected  waves  from  other  boundaries  of  the  metal  plate.  A  time 
of  about  7S-80|xs  passes  between  the  time  of  inq)act  and  the  time  of  crack  initiation.  In  this 
time  no  waves  traveling  in  the  PMMA  have  reached  the  notch  tip.  All  deformation  effects 
visible  in  the  interferograms  of  the  PMMA  side  are  from  leakage  of  energy  from  the  metal 
side  to  die  PMMA  through  waves  crossing  the  interface.  Essentially  throughout  all  tests, 
the  crack  is  driven  by  the  metal  side  of  the  bond.  This  will  become  even  more  obvious  in 
the  subsequent  analysis  of  the  high  speed  gas  gun  eiqieriments. 

After  the  crack  initiates  intense  stress  waves  are  emanated  from  die  crack  dp.  These 
waves  are  visible  in  figures  3  and  4  as  discrete  kinks  in  otherwise  smooth  fringes.  This 
observation  is  a  reliable  sign  of  the  existence  of  vary  strong  dynamic  effects  and  of  high 
velocities  of  crack  growth  that  ate  a  large  fraction  of  the  Rayleigh  wave  speed  of  PMMA 
(Cg*^).  It  should  be  noted  that  such  waves  are  not  usually  visible  after  crack  initiation  in 
homogeneous  PMMA  specimens  of  the  same  configuration  and  loading,  where  the 
observed  maximum  crack  tip  speeds  are  around  0.35c]^"*“.  However,  such  waves  do 
appear  in  optical  CXjS  or  caustic  patterns  of  growth  in  high  strength  materials  (see  Rosakis 
(1993)),  when  crack  tip  speeds  are  in  excess  of  0.4c™“*. 

The  actual  crack  tip  speed  and  acceleration  histories  corresponding  to  the 
PMMA/steel  experiment  shown  in  figure  4  are  plotted  in  figures  S(a,b).  The  crack  tip  speed 
increases  to  its  largest  value  in  a  very  short  time  (~15m.s).  This  is  also  evident  by  the 
presence  of  large  crack  tip  accelerations  at  the  beginning  of  crack  growth.  Consistent  with 
the  observation  of  stress  waves  during  crack  growth  is  the  measured  maximum  crack  tip 
speed.  This  is  seen  to  be  around  SOOm/s  or  0.85c]J*““.  This  observation  was  rather 
surprising  given  previous  experience  with  dynamic  crack  growth  in  homogeneous  PMMA 
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speciinens  of  the  same  configuration  and  loading  (maximum  speeds  of  ~0.3SCx^). 
Attributing  the  phenomentm  of  high  crack  dp  speeds  to  the  existence  of  a  weak  bond  is  not 
very  convincing  since  the  bond  calibration  technique  revealed  a  bond  fracture  toughness 
very  close  to  that  of  PMMA  (see  Tippur  and  Rosakis  (1991)).  In  fact  a  strong  bond  is 
necessary  in  ordn  to  drive  the  crack  tip  to  higher  speeds.  It  is  already  obvious  in  these 
experiments  that  energy  transfer  across  the  interface  fiom  the  metal  side  is  responsible  for 
crack  initiation  and  growth,  sinqrly  because  in  the  time  frames  involved  very  little  wave 
nootion  occurs  in  PMMA.  A  strong  bond  is  necessary  because  otherwise  less  energy  would 
be  transferred  and  the  crack  tip  would  not  be  able  to  propagate  as  fast.  This  fact  will 
become  even  more  obvious  in  the  higher  loading  rate  experiments  described  later.  In 
addition  we  believe  that  tiiese  effects  of  energy  transfer  are  mote  pronounced  because  of  the 
large  mismatch  of  elastic  properties  between  PMMA  and  steel  (or  aluminum). 

One  noteworthy  point  is  that  the  values  of  acceleration  quoted  in  figure  18(b)  are 
not  meant  to  be  exact  quantitative  results.  They  are  products  of  two  differentiations  of  the 
crack  length  history  and  are  susceptible  to  error.  Nevertheless  these  values  shed  light  into 
tire  orders  of  magnitude  of  acceleration  encountered  in  such  an  experiment 

Gas  gun  e^eriments 

Because  of  the  intrinsic  difficulties  associated  with  the  problem  of  dynamic 
interfacial  fracture,  few  theoretical  studies  have  bear  performed  (Goldshtein  (1967),  Brock 
and  Achenbach  (1973),  Atkinson  (1977),  BWillis  (1971,1973)).  Even  so,  it  should  be 
noted  that  there  is  a  disagreement  on  the  theoretically  predicted  value  of  the  terminal 
velocity  for  dynamically  propagating  interfacial  cracks.  Atkinson  (1977)  claimed  a  terminal 
velocity  equal  to  the  lower  of  the  Rayleigh  wave  speeds  of  the  two  constituents.  In  contrast 
Willis  (1973)  argued  that  the  terminal  velocity  was  a  little  greater  than  the  lower  of  the  two 
Rayleigh  wave  speeds.  More  recently  the  problem  of  dynamic  interfacial  crack  growth  has 
been  studied  by  Yang  et  al.  (1991).  They  provided  an  asymptotic  solution  for  the  stress 
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field  surrounding  die  tip  of  a  dynamically  propagating  Ismaterial  crack.  One  very  important 
conclusion  of  their  study  is  that  as  the  crack  dp  velocity  approaches  the  lower  of  the  two 
Rayleigh  wave  speeds  a  finite  amount  of  energy  has  to  be  transmitted  to  the  crack  tip  to 
maintain  extension  at  ^  with  non-zero  complex  stress  intensity  factor.  This  is  unlike  the 
homogeneous  case  where  an  infinite  amount  of  energy  has  to  be  transmitted  to  the  crack  tip 
to  maintain  extension  at  cg  if  the  stress  intensity  factor  is  non-zero  (Broberg  (1960)  and 
Freund  (1990)).  This  obviously  makes  it  energetically  impossible  for  a  crack  in  a 
homogeneous  solid  to  exceed  the  material’s  Rayleigh  wave  speed,  but  apparently  no  such 
restriction  «dsts  in  the  bunateiial  case. 

As  was  seen  in  figure  S,  experimental  measurements  from  the  drop  tower  tests 
resulted  in  crack  tip  speeds  close  to  (the  lower  Rayleigh  wave  speed).  In  fact 
towards  the  end  of  our  observation  window,  the  velocity  seemed  to  rise  further.  Since  no 
energetic  restrictions  exist  for  the  bimaterial  case,  we  believed  that  it  would  be  possible  to 
drive  the  interfacial  crack  to  higher  speeds  by  subjecting  the  specimen  to  different  loading 
histories  resulting  in  higher  rates  of  loading.  Use  of  a  high  speed  gas  gun  and  a  one  point 
bend  loading  configuration  were  made  to  produce  loading  rates  in  cracked  bimaterial  plates 
that  were  higher  than  those  encountered  in  drop  weight  tower  experiments.  The  gas  gun 
projectile  was  made  out  of  steel.  It  was  SOmm  in  diameter  and  75mm  in  length  and  the 
impact  velocity  used  was  2()m/s.  A  sequence  of  selected  interferograms  from  such  a  test  is 
shown  in  figure  6.  In  this  test  a  sharp  starter  crack  was  used.  The  resulting  crack  tip  speed 
and  acceleration  histories  for  this  particular  test  are  shown  in  figures  7(a,b).  It  is  clear  that 
the  Rayleigh  wave  speed  of  Plexiglas  (Cg*^)  has  been  exceeded  by  the  crack  tip.  In 
addition  the  shear  wave  speed  of  Plexiglas  has  also  been  exceeded.  After  the  shear  wave 
speed  of  PMMA  has  been  exceeded  the  crack  seems  to  decelerate  to  a  constant  speed 
around  1  lOOm/s.  This  feature  of  the  plateau  in  the  crack  tip  speed  history  was  also  seen  in 
several  other  gas  gun  tests.  The  reason  for  this  deceleration  and  constant  speed  is  not  yet 
clear  and  is  still  being  investigated.  In  any  case  it  is  a  v^y  repeatable  feature  of  the 
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experimental  measurements.  At  later  times,  though,  the  crack  tip  speed  increases  even 
fiirdier,  but  unfcntunately  the  growing  crack  exits  die  window  of  our  observation. 

When  a  specimen  containing  a  blunt  starter  notch  was  impacted  with  the  gas  gun 
projectile,  recorded  terminal  crack  tip  speeds  were  even  higher  than  those  seen  in  figure  7; 
in  some  cases  approaching  the  longitudinal  wave  speed  of  PMMA.  A  sequence  of 
interferograms  from  such  an  experiment  is  shown  in  figure  8.  The  corresponding  crack  tip 
speed  history  for  this  test  is  shown  in  figure  9.  Here  the  maximum  crack  tip  speed  is 
estimated  to  be  about  85%-90%  of  the  plane  stress  dilatadonal  wave  speed  of  PMMA. 
Visual  evidence  of  crack  propagation  in  the  intersonic  regime  for  PMMA  can  also  be  seen 
in  the  interferograms  in  figure  8.  (Recall  that  only  the  PMMA  side,  which  is  transparent,  is 
visible).  A  clear  change  of  the  nature  of  the  fringe  pattern  occurs  between  the  first  and  last 
frames  of  figure  8.  Initially,  in  the  first  one  or  two  frames,  the  fringes  are  smooth  and 
continuous.  Within  a  few  microseconds,  the  main  lobe  of  the  fringe  pattern  becomes 
squeezed  and  the  fringes  develop  distinct  kinks  along  certain  directions.  In  the  last 
interferogram  of  figure  8  the  kinks  have  merged  into  “shock  like”  formations.  It  is  not  clear 
at  this  time  exactly  what  these  kinks  represent  Since  the  crack  is  traveling  at  a  velocity  in 
the  intersonic  range  for  PMMA,  we  would  expect  shear,  but  not  dilatational,  shock  waves 
to  be  present  Unfortunately  because  the  optical  technique  of  CGS  measures  gradients  of 
direct  stress  components,  shear  effects  cannot  be  imaged.  It  is  believed  that  the  apparent 
“shock”  formations  in  the  later  interferograms  either  form  because  of  Ifead  waves,  that  are 
propagating  in  fiont  of  the  crack  tip,  or  are  three  dimensional  effects  associated  with  non- 
straight  propagation  of  the  crack  front  This  subject  is  still  under  investigation  and  will  be 
presented  in  Lambros  et  al.  (1994). 

These  observations  are  very  int^esting  because,  to  our  knowledge,  no  direct 
observations  of  unassisted  intersonic  or  supersonic  crack  propagation  have  ever  been  made 
in  homogeneous  materials,  even  though  a  large  number  of  theoretical  studies  exist  on  the 
subject  (e.g.,  the  contributions  of  Burridge  et  al.  (1979)  and  more  recently  of  0Broberg 
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(1985,  1989),  Bykovtsev  and  Kramarovskii  (1989),  Aleksandrov  and  Smetanin  (1990) 
and  Broberg  (1993)).  In  the  work  of  Broberg  (1985)  some  indirect  evidence  is  presented 
diat  shows  Aat  in  certain  eardiquake  sliding  events,  intersonic  crack  growth  velocities  may 
be  possible. 

In  experiments  involving  gas  gun  loading  the  observed  crack  initiation  time  was 
about  30iis  after  impact  This  is  much  less  titan  the  80p.s  of  the  previous  drop  weight  tower 
experimrats.  It  also  correlates  well  with  the  time  it  takes  a  conqtressive  dilatational  wave  to 
travel  the  width  of  tiie  specimen,  reflect  off  the  traction  free  bottom  surface  of  the  metal  side 
and  reach  the  tip  once  again.  This  reinforces  the  argument  that  the  whole  crack  growth 
process  in  tiiese  experiments  is  driven  by  the  metal  side  of  the  bond,  via  leakage  of  energy 
across  the  interface.  It  is  believed  that  it  is  an  excess  supply  of  this  energy  that  causes  the 
crack  tip  to  be  driven  at  speeds  that  are  intersonic  for  PMMA,  but  still  below  tite  Rayleigh 
wave  speed  of  steel. 

3.2  Data  reduction  and  results 

In  the  work  of  Tippur  and  Rosakis  (1991)  use  of  the  quasi-static  plane  stress 
asymptotic  stress  field  was  made  when  extracting  data  from  dynamic  interferograms. 
Obviously  this  is  not  a  correct  approach,  since  it  is  clear  that  large  dynamic  and  transient 
effects  occur  in  dynamic  interfacial  fracture.  However,  that  was  the  only  theoretical  result 
available  at  the  time.  Since  then  a  few  researchers  (Y ang  et  al.  (1991),  Wu  (1991),  BDoig 
(1992,1993))  have  determined  the  spatial  structure  of  the  near  tip  stress  field  around  a 
bhnaterial  crack,  propagating  dynamically  under  plane  stress  and  steady  state  conditions.  In 
addition,  a  higher  order  transient  analysis  for  a  dynamically  propagating  crack  has  been 
obtained  by  Liu  et  al.  (1993).  As  may  be  necessary,  one  can  use  either  the  asymptotic  or 
the  higher  order  transient  stress  fields  to  extract  the  relevant  fracture  parameters.  The 
procedures  have  been  described  in  Rosakis  et  al.  (1991)  and  Liu  et  al.  (1993).  Here  we 
shall  discuss  briefly  the  fitting  procedure  using  the  asymptotic  result  of  Yang  et  al.  (1991) 
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which  can  be  used  with  interferograms  exhibiting  some  region  of  dynamic  AT'^-dominance. 
(A  region  surrounding  the  crack  tip  in  which  the  stress  field  is  well  described  by  die  most 
singular  term  of  the  asymptotic  expansion  for  stress  is  called  a  AT'^-dominant  region). 

For  inteifacial  cracks  propagating  dynamically  under  steady  state  conditions  in 
bimaterial  specimens  Yang  et  al.  (1991)  observed  that  near  the  crack  tip  the  stress  field 
assumes  the  form. 


^=|<jg(d,v)+/i»i<  <T2(d.v) . 


(5) 


where  r,  dare  polar  coordinates  of  a  coordmate  system  translating  with  the  crack  tip  at 
speed  V  and  -  Kf  +  iK^  is  the  complex  dynamic  stress  intensity  factor.  Analytical 
expressions  for  and  are  given  in  Yang  et  al.  (1991).  Unlike  the  quasi-static  case, 

the  osdllatoiy  index,  e,  which  is  a  material  mismatch  parameter,  is  now  a  fimction  of  crack 
tip  speed  as  well  as  of  the  moduli  of  the  materials  of  the  bimaterial  combination  (i.e., 
e  =  ^v)). 

The  field  quantity  of  interest  in  analyzing  the  CCS  fiinge  patterns  for  PMMA  is 
chd{Oji+ 622)13X1  (see  equation  (4)).  By  differentiating  equation  (5)  with  respect  to  xi 

and  after  some  algebraic  manipulation,  we  have 

-(i  +  of  +  277a,  )cos^^  +  0  +  elnr^  (6) 

+2e(l +a^-  ^  -  e/n  r  j 

^  30 

-2e{l +al+  2Tia^)sin\^ + 0 + £//i  r 


where 
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(4a^,  -  (7 + oif^cosHea) 


fXu  = 


V/2 


7  — 


.  6i  =  uin~‘[{a^,)lx,\  r,=-Jx^,+a!xl, 


IC\i)  -  Kfit)+ iKf{t),  ^t)  =  B»-'(jir/(»)/A:;(»)). 


C/4  are  the  longitudinal  and  shear  wave  speeds  respectively  and  One  note  to  be 

made  here  is  tfiat  to  obtain  the  stress  field  in  the  lower  nuuerial  we  need  to  substitute  -n  for 
K  in  equation  (6).  In  the  present  study  though  we  only  obtain  interferograms  from  the 
deformation  of  material-1  (PMMA)  so  this  need  not  be  d(»e. 

Parameters  £  =  ^v),  ti  =  tI(v),  are  functions  of  crack  tip  spee  ^  id  material 
properties,  e  and  17  are  two  mianatch  parameters  that  charactoize  the  interface.  They  are 
derived  explicitly  in  Yang  et  al.  (1991) .  Parameter  £  is  given  by 


£(v)« 


with  2. 


(7) 


where  p  is  the  second  Dundurs’  parameter  characterizing  the  bimaterial  problem  (see 
Dundurs  (1969)).  Functions  Hap  depend  on  the  elastic  properties  of  the  two  constituents  of 
die  interface  and  the  crack  tip  velocity  v.  They  are  expliddy  stated  in  Yang  et  al.  (1991). 

The  variation  of  e  and  ij  with  crack  tip  speed  is  shown  in  figures  10  and  1 1 
respectively.  Both  diese  figures  are  for  a  PMMA/steel  bimatoial  combination.  The  plots  are 
for  crack  tip  speeds  below  the  Rayleigh  wave  ^eed  of  Plexiglas,  which  is  the  regime  in 
which  equation  (5)  is  valid.  The  parameters  are  such  that  £->£,,  77  ->  7  as  v-^0  and 
equation  (S)  reduces  the  quasi-static  result  of  Rice  (1988)  for  v=0.  Eg  is  the  quasi-static 
oscillatory  index. 
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In  the  above  relation  only  the  values  of  the  dynamic  coi!q>lex  stress  intrasity  factor 
and  die  crack  tip  ^leed  are  undetermined  1^  die  asyn^totic  analysis.  The  crack  t^  qieed  is 
an  external  parameter  that  must  be  furnished  by  separate  experimental  measurements.  In 
our  case  diese  measurements  are  made  by  recording  the  crack  loigth  variation  with  time. 
We  can  use  equation  (6)  to  extract  as  a  function  of  time  from  optical  interferograms  by  a 
performing  a  least  squares  fitting  procedure.  Irom  the  above  discussion  it  becomes 
obvious  that  extraction  of  parametos  like  is  possible  provided  diat  eiqierimaital  data  are 
gathered  from  a  region  near  the  moving  crack  tip  characterized  by  die  structure  presented  in 
equations  (S)  and  (6),  i.e.,  that  there  exists  a  re^on  of  AT'^-dominance.  However  in  the  case 
that  no  region  in  die  specimen  is  characterized  by  such  a  AT'^-field,  the  use  of  a  hi^ier  order 
analysis  is  mandated  as  Ascribed  in  liu  et  al.  (1993).  The  least  squares  fitting  procedure 
utilized  in  this  study  is  very  similar  to  that  used  in  Mason  et  al.  (1992)  or  Rosakis  (1993) 
and  in  the  interest  of  brevity  will  not  be  described  here.  One  point  to  note  thou^  is  that 
when  analyang  optical  data,  care  must  be  takm  to  use  only  data  points  r/,  6Sf  from  a  region 
of  the  specimen  that  does  indeed  qiproximate  two  dimensional  plane  stress  conditions,  if 
we  wish  to  use  equations  (4)  and  (6)  to  interinet  the  experimental  results.  The  region  in 
which  such  conditions  are  valid  in  a  three  point  bend  bimatraial  specimen  was  determined 
by  the  finite  element  numerical  analysis  of  Lee  and  Rosakis  (1993).  Selection  of  data  points 
was  made  on  the  basis  of  the  results  of  this  analysis. 

The  variation  of  stress  intensity  factors  Kf  and  AT/  PMMA/steel  drop 

weight  tower  test  presented  in  figure  4,  are  shown  in  figures  12(a,b).  The  energy  release 
rate  G  and  phase  angle  <f>  for  the  same  test  are  shown  in  figures  13(a,b).  G  for  the  case  of 
dynamic  intorfadal  crack  growth  is  given  by  (Y ang  et  al.  (1991)), 

G  =  /-|Arf, 


(8) 
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Mdiere  F  is  a  function  of  velocity  and  biniaterial  system  material  properties  shown  e}q)licitly 
in  Yang  et  al.  (1991).  Note  titat  F  is  finite  at  the  Rayleigh  wave  speed  of  Plexiglas.  The 
phase  angle  in  figure  13(b)  has  been  calculated  fiom 


f 

V 


(9) 


This  implies  that  a  length  scale  of  Im  has  bem  used  in  this  particular  plot  The  phase  angle 
ranges  fitMn  -nto  nto  cover  all  possible  sign  combinations  of  Kf  and  ATj . 

On  all  plots  in  figures  12  and  13  t=Ofis  corte^nds  to  the  estimated  time  of  crack 
initiation.  Pre-initiation  values  of  Kf  and  Kf  were  detomined  by  fitting  the  quasi-static 
stress  field  of  Rice  (1988)  to  interferograms  before  t=0fu.  Unfortunately  determining  the 
time  of  initiation  in  each  test  is  not  an  absolute  science.  It  was  observed  tiiat  when  fitting  a 
K^-Geid  to  pictures  a  few  microseconds  after  r^O/is  the  fit  was  good  only  on  the  outer  most 
fringes.  This  is  contrary  to  what  would  be  expected  at  first  glance,  on  the  basis  of 
as3rmptotics.  An  exanq>le  of  such  a  case  is  shown  in  figure  14  for  the  interferogram 
corresponding  to  t=2fis  for  the  experiment  shown  in  figure  4.  The  diamond  points  are  data 
points  digitized  frran  the  interferogram,  and  the  solid  line  is  the  result  of  the  fit  of  equation 
(6)  to  those  same  points.  It  can  be  seen  that  the  three  outer  fringes  are  in  very  good 
agreement  with  the  ilT^-domiruuit  fit,  while  the  inside  fringes  are  not  Gradually  as  time 
progressed  (2-lOps)  the  quality  of  the  fit  deteriorated  from  the  inside  out  until  all  K^- 
dominance  was  lost  The  gradual  collapse  of  the  ilT^-field  between  pictures  agreed 
reasonably  well  timewise  with  a  dilatational  wave  being  emitted  from  the  crack  tip  and 
traveling  radially  outwards.  Such  waves  would  indeed  start  being  emitted  when  the  crack 
tip  initiated.  They  would  travel  outwards  and  destroy  the  IT^-field  that  was  established  in 
tile  pre-initiation  stages  of  the  eiqierimait  It  was  possible  to  determine  tiie  time  of  initiation 
by  choosing  the  interferogram  immediately  before  the  first  sign  that  the  innermost  firinge 
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did  not  conform  to  a  X'^-dominant  field,  Le.,  just  before  the  emission  of  a  wave  from  die 
crack  tip.  Even  so  there  is  some  uncertainty  in  the  exact  time  of  crack  initiation,  based  on 
when  the  crack  is  visually  seen  to  propagate.  This  uncertainty  is  represented  in  figure  12 
and  13  as  being  between  die  two  dotted  vertical  lines. 

It  can  be  seen  in  figure  12(a,b)  that  both  Kf  and  K2  increase  to  a  maximum  value 
and  dien  sharply  drop,  with  Kf  attaining  negative  values  towards  the  end  of  the  test  The 
energy  release  rate  (figure  13(a))  also  increases  to  a  peak  a  litde  after  r=0/£s  and  then  drops 
to  a  value  close  to  zero  implying  that  the  process  of  dynamic  crack  propagation  along 
interfaces  is  highly  unstable.  This  result  is  corroborated  by  the  steqi  rise  of  the  crack  tip 
speed  in  very  short  times.  The  continuing  increase  of  G  for  a  short  time  after  crack 
initiation  is  lelaied  to  die  fact  diat  a  material  particle  away  from  die  crack  tip  is  not  aware  of 
initiation  until  stress  waves  entitled  from  the  tip  reach  it  This  effect  was  illustrated  in  figure 
14.  It  has  also  been  observed  by  other  researchers.  For  exanqile,  Aoki  and  Kimura  (1993) 
in  a  numerical  investigation  using  synthetic  caustic  patterns  at  the  tip  of  a  crack,  observed 
an  overshoot  of  the  measured  G  just  after  initiation.  On  the  other  hand  the  phase  angle  0 
increases  throughout  the  test  In  equation  (6)  it  was  seen  that  in  the  theoretically  predicted 
stress  field  from  the  Yang  et  aL  (1991)  AT'^-dominant  analysis  velocity  v  and  phase  angle  0 
have  conqilimentary  effects  since  they  are  both  inside  the  argument  of  die  cosine  function, 
(v  enters  the  cosine  through  the  velocity  dependence  of  e).  Since  both  v  and  0  are  seen  to 
vary  considerably  throughout  crack  propagation,  but  the  size  and  orientation  of  the 
observed  fringe  pattern  around  the  tip  seems  unchanged  (see  figures  3  and  4),  the  two 
quantities  have  to  be  linked  in  some  way  as  to  keep  the  stress  gradients  resulting  from 
equation  (6)  constant  This  fact  will  be  exploited  further  in  the  next  section  to  propose  a 
possible  dynamic  crack  growdi  criterion  for  interfacial  fracture. 

Because  of  the  presence  of  oscillatory  tmns  in  the  stress  field  shown  in  equation 
(6)  all  inconsistoicies  associated  with  die  the<^  of  static  linear  elastic  interfacial  mechanics 
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are  still  piesent  Of  paramount  importance  is  die  inconsistency  of  predicted  inteipenetratitm 
of  the  crack  faces.  The  «ctait  of  the  resulting  contact  zone  is  given 

r^-e - c — (10) 

To  accept  the  measured  values  of  G  and  ^  as  accurate  we  must  make  sure  tiiat  the  extent  of 
predicted  contact  is  much  smaller  than  the  smallest  characteristic  length  of  the  problem  and 
is  embedded  in  die  three  dimensional  region  of  deformation.  (Dnly  in  that  case  can  die  K^- 
field  surrounding  the  crack  t4>  control  die  firacture  process.  From  figure  10,  which  shows 
the  variation  of  e  widi  crack  tip  qieed  for  a  PMMA/steel  system,  it  can  be  seen  that  as 
V  ->  e  400.  As  V,  and  consequently  £,  gets  larger  it  would  seem  that  the 
predicted  contact  zone  would  become  unaccqitably  large.  When  calculating  though 
from  experimental  measurements  it  was  found  that  die  variation  of  (jf  suppressed  the 
increase  of  the  contact  zone  size  keeping  it  at  very  low  values  for  all  but  the  highest  crack 
tip  velocities.  An  idea  of  the  size  of  the  contact  zone  predicted  for  the  one  point  bend 
PMMA/steel  test  of  figure  4,  is  givoi  in  figure  15.  The  solid  lines  are  contours  of  constant 
contact  zone  as  derived  from  equation  (10)  for  a  range  of  v  and  0.  The  diamor^  points  are 
the  calculated  experimental  results.  Clearly  values  of  the  contact  zone  are  below  10^. 
This  is  well  within  the  zone  of  three  dimensional  deformation  surrounding  the  crack  tip 
whose  maximum  extent  is  artmnd  2mm. 


4.  A  criterion  for  dynamic  interfadal  crack  growth 

A  sequence  of  dynamic  interferograms  obtained  from  the  tip  of  a  crack  lying  along  a 
PMMA/aluminum  interface  was  shown  in  figure  3.  Pre-  and  post-  initiation  pictures  were 
^own.  As  was  pointed  earlier,  before  crack  initiation  {KOfis)  the  fringe  pattern  around  the 
tip  changes  both  in  size  and  orientation,  denoting  a  change  in  the  energy  release  rate  G  and 


-18- 


phase  angle  0.  A  variation  was  indeed  seen  in  the  measured  values  of  G  and  ^  at  times 
before  for  die  I^fMA/steel  drop  weight  tower  test  seen  in  figure  4.  In  contrast,  aftor 
initiation,  the  fringe  pattern  surrounding  the  crack  does  not  change  significandy  with  time. 
The  implicadon  of  this  visual  observation  is  that  some  fundamental  physical  quantity,  such 
as  stress  or  crack  opening  diqdacement,  must  ronain  constant  throughout  die  crack  growth 
phase.  It  has  already  been  established  that  v,  ^  and  G  (or  all  vary  considerably  during 
crack  growdi.  Therefore,  if  we  wish  to  mathemadcally  model  die  fact  that  the  fringe  pattern 
during  crack  growth  does  not  change,  we  must  employ  some  combination  of  these  three 
quantities  diat  remains  constant 

It  was  seen  from  figure  15  that  die  variations  of  v  and  ^  were  such  as  to  keep  the 
size  of  the  predicted  contact  zone  small  throughout  each  test  This  suggests  that  there  is  a 
sqiarate  link  between  v  and  ^during  crack  growdi.  In  figure  IS  it  is  clearly  visible  that  die 
diamond  points  (corresponding  to  experimental  results)  curve  to  the  right  denoting  an 
increase  in  both  v  and  ^  as  the  crack  grows.  We  intend  to  investigate  whethor  the  particular 
padi  of  V  and  ^  in  each  experiment  is  such  as  to  maintain  some  fundamental  quantity 
constant  A  candidate  for  such  a  quantit>'  would  be  something  that  does  not  involve  \K^. 
We  will  look  into  the  ratio  of  the  opening  (^)  to  shearing  {5i)  displacements  behind  the 
crack  tip  as  one  possible  candidate,  ^th  referaice  to  the  geometry  of  figure  2  Sj  and  82  are 
d^nedas 


8Sr)-u^{r,e  =  +it)-u^{r,e--n),  a  =  1,2,  (11) 

where  Uoi(r,0)  is  the  displacement  field  surrounding  the  crack  tip.  Figure  16  shows  a 
geometrical  representation  of  8j  and  82  behind  the  crack  tip.  A  positive  82  means  an 
opening  displacement  and  a  positive  81  represents  sliding  of  material- 1  along  the  positive 
X2  axis.  From  the  asymptotic  analysis  of  Yang  et  al.  (1991),  the  crack  face  displacements 
are  given  by. 
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2r 

L 

^  1. 

/  +  4£^ 

siii(^  +elnr-  tan~'(2e)) , 


and 


(12) 


Quantities  e  and  rj  have  been  previously  defined  and  the  function  H22,  which  is  a  function 
of  mattrial  properties  and  crack  tip  speed,  is  given  in  Yang  et  al.  (1991).  Each  of  the  two 
displaconents  in  equation  (12)  depends  <xi  \K\  but  their  ratio  is  (»ly  a  function  of  v  and  (ft. 
It  is  given  by. 


J^(r)  =  ^tan{(p +elnr-  /a«~'(2£)) 


(13) 


If  we  assume  that  throughout  propa^on  dus  ratio  ranains  constant,  say  at  a  value 
Cj,  at  a  fixed  distance,  say  o,  behind  the  crack  t^,  then  equation  (13)  defines  a  single  line 
in  the  V,  ^plane.  The  equation  of  this  line  is  givoi  by, 

S  1 

~Cj=  ~~tan{(p +elna-  tarr^{2e)), 

or  solving  for  0, 

0  =  tarT^TiCj) + tan~\2e)  -elna.  (14) 


The  dependence  of  ^  on  v  comes  in  through  the  velocity  dependence  of  the  oscillatory 
index  e  and  the  traction  resolution  factor  t}.  (See  figures  10  and  1 1  for  the  variation  of  e  and 
rj  with  V  req)ectively). 
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In  figure  17,  which  is  a  plot  in  the  (p,  v)  plane,  the  experimental  results  of  die 
measured  crack  tip  speed  v  and  calculated  phase  angle  ^  are  shown  as  diamonds.  These 
results  are  fnrni  the  PMMA/steel  one  point  bend  experiment  conducted  under  drop  weight 
tower  impact  at  4m/s  which  was  illustrated  in  figure  4.  The  solid  line  in  figure  17 
represents  the  curve  given  by  equation  (14)  that  best  fits  the  experimental  values.  It  was 
decided  to  evaluate  equation  (14)  at  a  distance  a-2mm,  which  is  approximately  the 
innermost  radius  of  the  region  of  JlT^-dominance  observed  in  most  interferograms.  This 
means  that  it  is  also  the  point  closest  to  the  crack  tip  where  equation  (14)  is  valid.  The 
confuted  ratio  SitSi  for  this  fit  is  Cj=-03.  Bom  figure  17  it  seems  diat  indeed  throughout 
crack  propagation  the  quantity  81/82  at  2mm  behind  the  moving  crack  tip  remains  a 
constant  The  amount  of  crack  face  opening  is  3.3  times  the  amount  of  crack  face  sliding, 
dius  indicating  growdi  at  a  primarily  qiening  mode  at  this  particular  distance  behind  the  tip. 
The  choice  of  a  as  being  2mm  is  not  cracial  for  tiiis  result  At  a  difiinait  distance  behind  the 
tip,  the  ratio  of  shearing  to  opening  displacements  would  again  be  constant  Because  of  the 
nature  of  the  intofacial  field,  however,  the  numerical  value  of  this  constant  would  not  be 
the  same.  But  the  results  presented  in  this  weak  would  still  hold. 

The  condition  of  a  constant  81/82  behind  the  crack  tip,  equation  (14),  cannot 
constitute  a  dynamic  crack  growth  criterion  on  its  own  because  it  does  not  predict  the 
magnitude  of  the  stresses  surrounding  the  propagating  crack.  It  may  form  part  of  a 
propagation  criterion  that  relates  the  phase  angle  with  the  crack  tip  speed.  The  size  of  the 
stress  field,  and  the  CXjS  fringe  pattern,  surrounding  the  crack  tip  is  primarily  dependent 
upon  \K^\,  the  magnitude  of  the  complex  stress  intensity  factor.  The  second  part  of  a 
fracture  criterion  would  thus  have  to  involve  a  link  between  \K<^\  (or  G),  v  and  (jf.  Consider 
a  constant  crack  opening  displacement,  say  of  value  C2,  at  the  same  distance  a  behind  the 
crack  tip  as  before.  From  equation  (12b)  this  can  be  expressed  as 
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52^  = 


cosk^xe) 


cos{^ + elna  -  tarr^{2£))  =  Cj, 


or  solving  for 


‘>ll-¥4e^  cosh{n£)  jlT  1 

Hj2  1 2a  cos{(p + elna  -  tan~\2e)) 


(15) 


The  above  is  a  relation  between  [K^\,  v  and  0  with  C2  and  a  as  parameters.  For  specific 
values  of  C2  and  a  equation  (IS)  represents  a  surface  in  the  UT^t,  v  and  0  space.  An  example 
of  such  a  surface  in  the  G,  v,  ^  space  (by  using  equation  (8))  is  plotted  in  figure  18  with 
a=2mm  and  C2=lNlm^  (arbitrary). 

We  can  now  combine  equations  (14)  and  (IS)  to  propose  a  conq>lete  criterion  for 
dynamic  crack  growth  along  bunaterial  interfaces.  Substituting  in  equation  (IS)  for  0  from 
equation  (14)  we  obtain  a  relation  between  l&f  and  v  only  which  incorptnates  the  fact  that 
S]/d2  is  constant  The  combination  of  (14)  and  (IS)  yields, 


|r'|=c, 


cosMjne)  _ [ _ 

H22  i2a  cos(fan~^(CjTi)y 


(16) 


The  above  depends  on  constants  C/  and  C2.  Dependence  on  parameter  a  comes  in  only 
through  its  necessity  in  the  choice  of  a  constant  Cj  to  fit  the  experimental  data  using 
equation  (14),  as  was  done  in  figure  17.  The  value  of  Ci  is  then  fixed  after  performing  this 
fit  For  the  results  shown  in  figure  17  it  was  seen  that  Ci=-03.  Unfortunately  constant  C2, 
which  essentially  is  the  opening  strength  of  the  bond,  carmot  be  easily  determined.  It  can  be 
eliminated  though  by  normalizing  equation  (16)  with  some  reference  value.  If  we  assume 
that  one  of  the  experimentally  obtained  data  points,  with  parameters  vo,  eo,  t\o,  ^22  [^0 1 
and  Go,  follows  equation  (16)  then  we  can  eliminate  C2  by  conq)uting  |Ar''|/|Ar0 1  as 
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^(j + 4^)cosh{nE)H22  cos(tan~*{CjiJo)) 
^(l + 4^)cosh{iteo)^22  cos(tan~\C,Ti)) 


(17) 


The  above  result  can  be  converted  to  a  ratio  of  energy  release  rates  G/Gq  u^g  equation 
(8).  Comparison  with  all  other  experimental  points  is  then  possible  by  plotting  the 
equivalent  of  equation  (17)  in  die  G,  v  plane. 

For  tire  PMMA/steel  drop  weight  tower  test  used  previously,  a  comparison  of  the 
experimental  data  with  equation  (17)  is  shown  in  figure  19.  The  normalizing  point  used  is 
die  mre  corresponding  to  v^lSOmls.  This  point  is  forced  to  agree  with  die  computed  curve 
(solid  line  in  figure  19)  since  it  has  been  used  to  eliminate  the  constant  C2.  Nevertheless, 
the  remaining  points  are  not  forced  to  follow  die  curve.  It  can  be  seen  though  that 
agreement  is  very  good.  A  notable  issue  hoe  is  that  die  line  in  figure  19  is  /mr  a  fit  to  die 
expaimental  values.  Only  Ci  is  fitted  to  the  experimental  values  using  equation  (14)  as 
described  above.  Once  this  is  done  the  curve  given  by  equation  (17)  is  completely 
qrecified.  This  implies  that  indeed  not  only  is  the  ratio  dilSz  constant  behind  the  crack  tip, 
but  82  itself  is  constant  These  two  facts  therefore  constitute  the  proposed  dynamic  fracture 
critetitm  for  interfaces.  An  inteifacial  crack  that  has  beoi  impaa  loaded  will  initiate  dynamic 
propagatitm  at  certain  critical  values  of  Gc  and  as  is  sug^sted  by  the  work  of  Liechti  and 
Qiai  (1991;  Liechti  and  Qiai  (1992).  Subsequoitly  propagation  will  be.such  as  to  maintain 
Si! 82  and^  constantatacriticaldistancebehindthemovingcracktip.  Alternatively,  the 
prqiosed  criterion  can  be  viewed  as  critical  opening  (82  )  and  shearing  (81)  displacements 
at  some  distance  a  behind  die  moving  crack  tip.  This  is  equivalent  to  requiring  a  self-similar 
crack  profile  throughout  crack  growdi. 

A  possible  path  in  G,  v,  ^  space  followed  during  interfacial  crack  growth  that  obeys 
the  above  proposed  criterion  is  seen  in  figure  18.  As  was  mentioned  above,  this  surface 
represents  the  result  of  equation  (IS).  For  v=Omfs  it  reduces  to  the  usual  “U"  shaped  curve 
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that  has  been  observed  in  quasi-stadc  experiments  by  Liechd  and  Chai  (1991).  This  serves 
as  an  initiation  condition  for  the  case  of  impact  loading.  The  currently  proposed  criterion 
therrfore  will  act  as  follows : 

1. -  When  the  stress  field  surrounding  the  stationary  crack  tip  reaches  a  critical  value 
of  Gc  and  ^  (or  equivalently  Sj  and  then  crack  initiation  occurs. 

2. -  Subsequently  the  dynamically  growing  crack  will  always  propagate  in  such  a 
way  as  to  keep  Sj(t)  and  ^2(0  behind  it  constant  at  their  initiation  values.  This  is 
equivalent  to  maintaining  a  constant  crack  profile.  Such  a  path  is  shown  on  the  surface  of 
figure  1 8  as  die  line  on  the  surface  curving  to  the  tight  The  projection  of  this  line  on  the  v, 
^  plane  will  be  the  contour  of  constant  Sj /Sj  seen  in  figure  17. 

Finally,  a  very  inqwrtant  point  to  note  is  that  although  the  proposed  criterion  was 
derived  using  the  results  of  the  complex  stress  intensity  factor  history,  the  criterion  itself 
may  be  general  enough  not  to  require  i^^^-dominance  in  its  application.  The  criterion  of  a 
critical  crack  tip  opening  displacement,  or  opening  angle,  has  been  used  extensively  in  the 
past  for  the  fracture  of  homogeneous  solids.  It  has  been  seen  to  be  equivalent  with  initiation 
criteria  based  on  Kj  or  7/  and  growth  criteria  (quasi-static  or  dynamic)  based  on  Kj  in 
situations  where  liT-dominance  was  assured.  However,  a  criterion  based  on  crack  tip 
opening  displacement  has  also  been  found  to  control  fracture  in  instances  when  Kj  orJj  are 
invalid  or  cannot  be  defined,  such  as  growth  in  elastic-plastic  materials.  We  therefore 
strongly  believe  that  the  proposed  criterion  for  dynamic  interfacial  crack  growth  based  on 
both  opening  and  shearing  displacements  behind  the  crack  tip  (crack  profile)  will  be  valid  in 
a  wide  array  of  cases  as  well;  even  if  a  /ir<^-dominant  region  cannot  be  found. 

In  general,  though  the  particular  numerical  values  obtained  in  this  study  should  be 
critically  accented,  the  trends  indicated  seem  well  founded  and  logical.  Clearly  more  tests 
are  required,  using  a  wider  array  of  loading  geometries  and  loading  rates,  to  obtain  a 
complete  understanding  of  the  phenomena  involved  in  dynamic  interfacial  fiacture.  This  is 
currently  underway  and  the  results  will  be  reported  in  a  future  publication  dealing 
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exclusively  with  intofacial  crack  growth  criteria.  Nevertheless,  we  believe  that  a  dynamic 
prcqMigation  criterion  involving  critical  displacemoits  behind  the  crack  tip  is  both  reasonable 
and  accurate,  and  should  be  investigated  further.  For  a  more  extensive  discussion  of  the 
proposed  fracture  criterion  see  Lambros  and  Rosakis  (1994). 

The  concept  of  a  crack  growth  criterion  for  dynamic  interfacial  fracture  has  also 
been  investigated  numerically  by  Lo  et  al.  (1993).  More  recent  results  of  their  work,  that 
appear  in  the  same  special  issue  of  this  journal,  discuss  a  comparison  of  the  current 
experimental  results  to  a  numerical  simulation  of  our  experiments  (see  Nakamura  et  al. 
(1994)). 


5.  Conclusion 

Very  large  dynamic  effects  were  observed  in  all  dynamic  bimaterial  tests.  In  the 
case  of  the  high  speed  gas  gun  experiments,  a  terminal  crack  propagation  velocity  around 
90%  of  the  plane  stress  dilatational  wave  speed  of  PMMA  was  observed.  Because  of  the 
very  short  duration  of  the  experiments  in  this  study,  it  may  be  concluded  that  the  whole 
process  of  interfacial  crack  initiation  and  growth  in  these  tests  is  driven  by  energy  ’leaking” 
from  the  iTKtal  side  to  die  PMMA  side  of  the  bond.  It  is  believed  that  this  fact,  enhanced  by 
the  large  elastic  property  mismatch  of  the  constituents  of  the  bimaterial  system,  is  what 
causes  the  high  dynamic  nature  of  these  tests.  In  addition  to  dynamic  effects,  very  severe 
transient  effects  have  been  seen  to  occur,  especially  during  the  early  stages  of  crack 
growth.  Evidence  of  transience  are  the  very  large  accelerations  recorded  (maximum  values 
of  the  order  of  lO^g  where  g  is  the  acceleration  of  gravity).  The  phenomenon  of  the 
existence  of  highly  transient  deformations  was  corroborated  by  direct  comparison  of  the 
experimental  results  to  a  transient  elastodynamic  analysis  in  the  work  of  liu  et  al.  (1993). 

It  was  experimentally  observed  that  the  crack  profile  (opening  and  shearing 
displacements  at  a  given  distance  behind  the  propagating  crack  tip)  remained  constant 
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throughout  propagation.  This  observation  formed  the  basis  of  the  proposed  dynamic 
fracture  criterion.  The  use  of  the  notion  of  an  invariant  crack  face  profile  as  a  fracture 
criterion  is  both  physically  sensible  and  consistent  widi  fracture  models  prqwsed  for  quasi¬ 
static  and  dynamic  fracture  of  homogenecMis  materials.  Even  tlK>ugh  in  this  particular  study 
use  of  die  coirqilex  dynamic  stress  intensity  factor  was  made  in  the  observation  of  a  critical 
state  in  crack  face  displacements,  the  proposed  fracture  criterion  itself  is  of  a  more  general 
nature.  A  criterion  based  <hi  crack  face  profile  may  still  be  tqiplicable  to  situations  where 
large  plasticity  and/or  large  scale  diree  dimensional  effects  are  dominant  or  where  intersonic 
and  supersonic  crack  growth  is  occurring  (e.g.,  in  a  situation  where  there  is  no  K*^- 
dominance). 
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^peed  Camera 


Figure  1:  Schematic  of  CGS  set-up  in  transmission. 


Figure  2;  Geometry  of  interfadal  crack  problem. 
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DYNAMIC  CRACK  GROWTH 


t=0/is  t=14/is  t=2S//s 


Figure  3  :  Selected  sequence  of  CGS  interferograms  from  a  PMMA/aluminum 
three  point  bead  drop  weight  tower  test  (from  Tippur  and  Rosakis  (1991)). 


PMMA/steel  drop  weight  tower  test. 
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Figure  6:  Selected  sequence  of  CGS  interferograms  obtained  from  a  one 
p<dnt  bend  air  gun  test  on  a  PMMA/steel  specimen. 


Figure  8:  Selected  sequence  of  CGS  interferograms  of  a  growing  crack  in  a 
one  point  bend  PMMA/steel  air  gun  experiment  (blunt  starter  notch). 
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t  Cus) 


Figure  9:  Velocity  time  history  for  the  experimoit  shown  in  figure  8. 
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Figure  10:  Variati<Mi  of  oscillatory  Index  with  velocity  for  a  PMMA/steel 
bimaterial  system  (plane  stress  conditions). 


Figure  11:  Variation  of  traction  resolution  factor  with  velocity  for  a 
PMMA/steel  bimaterial  system  (plane  stress  conditions). 
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Figure  13(a,b):  Time  history  of  energy  release  rate  and  phase  angle  (L=lm) 

for  the  experiment  shown  in  figure  4. 
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Figure  15:  Magnitude  (rf  predicted  contact  wot  in  PMMA/steel  one  point 
bend  drop  weight  tower  experiment  shown  in  figure  4. 


interface 

Figure  16:  Schematic  definition  of  opening  and  shearing  displacements 

behind  the  interfacial  crack  tip. 


Figure  18:  (G,  v,  f)  surface  for  the  quantity  =  const.  Crack  path 
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Figure  19:  Comparison  of  normalized  energy  release  rate  from  the 
PMMA/steel  drop  weight  tower  experimoit  with  the  prediction  of  constant 
crack  profile  growth  (equation  (17)). 
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PROPERTY 

FMMA 

6061-T6  Aluminum 

AISI 4340  Steel 

E(GPa) 

3.24 

80 

208 

V 

0.35 

0.33 

0.3 

2080 

6600 

5970 

IIIIIIIIIHSRIliSIIIH 

1000 

3330 

3190 

IliRR^HI 

935 

3100 

2950 

p  (Kgi/niP) 

1190 

2710 

7830 

=  0.mi  e‘‘'^  =  0.1037  e"*  =  0.1073 

Table  1:  Mechanical  properties  of  interface  constituents. 
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